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Abstract 

We propose a nonperturbative scheme for the calculation of thermal damping-rates using 
exact renormalization group (RG)-equations. Special emphasis is put on the thermal RG 
where first results for the rate were given in Q. We point out that in order to obtain 
a complete result that also reproduces the known perturbative behaviour one has to take 
into account effects that were neglected in [|J. We propose a well-defined way of doing 
the calculations that reproduces perturbation theory in lowest order but goes considerably 
beyond perturbative results and should be applicable also at second order phase-transitions. 
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Perturbation theory at finite temperatures is often invalidated by bad infrared behaviour. 
These problems can be solved by different resummation schemes. One of the most powerful 
approaches to resummation is the renormalization group of Wilson and others 0. This involves 
introducing an external scale and deriving functional differential equations for the dependence of 
generating functionals on this scale. The right hand side of such an RG-equation can formally 
be interpreted as a one-loop expression in the sense that it is of order h compared to the left 
hand side. Nevertheless the Wilsonian RG constitutes a nonperturbative method and the RG- 
equations are exact functional relations. The approach is well known to correctly reproduce the 
infrared-behaviour of theories even at second-order phase-transitions. 

Even though the RG-equations define a nonperturbative approach, it is possible to reproduce 
perturbation theory. This is straightforward for one-loop results, but becomes rather tedious 
if one goes to higher orders. Quantities which are well described perturbatively for a range of 
temperatures but for which perturbation theory fails e.g. at a phase-transition constitute a major 
challenge for RG-equations. This holds in particular for quantities where the one-loop contribution 
vanishes. The thermal damping-rate in <£> 4 -theory is an example for a quantity where the lowest 
order perturbative contribution is two-loop. If one aims at a reliable calculation of such quantities 
at all temperatures, one has to reproduce two-loop perturbation theory as the leading term and 
apply sensible nonperturbative resummations close to a possible phase-transition. 

A recent formulation of Wilsonian RG-equations in the Schwinger-Keldysh (CTP-) formalism 
of real-time thermal field theory is particularity suited for calculations of nonstatic thermal Green- 
functions ||. It makes use of the fact that the propagators in the CTP-formalism separate into 
the usual zero-temperature and a finite-T part. The finite-T parts only contribute on-shell and 
depend only on the three-momenta. The cutoff modifies the thermal part and one introduces 
cut-off propagators of the form 




°^ = [ -A- ) + < A " A * 

+(A-A')6(|S|,A)iV(|y) 

where 0(|fc|,A) is a possibly smeared out step function. We will in the following use a sharp 
cutoff, i.e. set 0(|fc|, A) = 0(|A;| — A). Thus for finite A, the propagation of thermal modes with 
three-momemtum small compared to A (soft modes) is supressed, while the hard thermal modes 
are unmodified. Inserting this propagator into the usual expression for the generating functional 
Z[J] one obtains a A-dependent functional Z^. Introducing the modified Legendre-transform 

r A [$] = -ilnZ A [J] - J • $ - 1$ • (D^- 1 • $ (2) 

one readily obtains for the scale-dependence of Ta the exact functional differential equation || 

A^I = Itva^ (V + WY 1 (3) 
OA 2 OA \ A 5<5>d$ / v ; 
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This is the thermal renormalization group equation (TRG) and will be the starting point of our 
discussion. 

Being a differential equation, (Q) of course has to be supplemented by boundary conditions. As 
discussed in ||, in the limit A — > oo, the effective action T\ is trivially obtained from the (T = In- 
effective action of the theory. Also, in the limit A — ► the full finite temperature CTP-effective 
action is obtained. 

RG-equations for Green-functions are simply obtained from (|3]) by taking functional derivatives 
with respect to $. For the damping-rate, we are interested in the imaginary part of the two-point 
function. The flow-equation for the two-point function reads in terms of the thermal fields tpi with 
i = 1, 2 (y? is some background configuration) 



A 9 S2T > 



dA SifiSifj 



7=0 



-TrK 

2 5ipi6<pj5$6$ 



•Pa 



2 5(pi8<f>8Q 5(pj5§5& 



+ (*<-> j) (4) 



¥0 

where K\ = —iVA 9 ^ V is the kernel of the TRG and T> = (j)^ 1 + is the full two-point 

function 0]. Even though in this representation the flow-equation appears like a plain one- loop 
equation, due to the fact that all n-point functions involved are full n-point functions it is actually 
an exact, nonperturbative expression. 

Let us then assume that we want to calculate the imaginary part of the thermal two-point 
function in a scalar theory with unbroken Z 2 -symmetry. Perturbatively in this case the lowest 
contribution to the imaginary part of the self-energy occurs on two-loop level. As out in in 
such a situation the contribution to the imaginary part of (|]) has to come from an imaginary 
part of the full four-point function (by the ^-symmetry, the second contribution in (|4]) vanishes 
identically for all A). We thus have to solve a coupled system of at least two flow-equations, (f|) 
and a corresponding equation for the imaginary part of the four-point function. 

First however let us make some remarks concerning the thermal indice which appear in real- 
time formulations of thermal field-theory [|J . The physical fields always carry a thermal index 1 
and one is in general interested in the calculation of Green-functions for those fields. Since the 
propagators are nondiagonal in the thermal indice, in the calculation of higher loop contributions 
one has to include vertice of the 2-fields nevertheless - neglecting those contributions yields singular 
expressions. If we use the TRG to compute Green-functions, we have to use full vertice according 
to (^). In this case, one has to allow for vertice with mixed thermal structure. 

Furthermore, what enters the physical self-energy (and thus the damping-rate) is not the 
Green- function S S J^ 1 , but rather a specific combination of the 0, known as fl(p) in the 

literature. Assuming a Schwinger-Dyson equation for the full propagator it is easy to show that 
for the real- and imaginary parts of IT the following relations hold: 

Tzu(p) = n 

III l( 62t 62t \ 

° \8<pi(p)5ipi(-p) Scpi(p)5cp 2 {-p) J 



3 



(e(po) is the sign function). XTI can also be obtained as 

i 5 2 r 

Tn( v ) = _ j — (q) 

{P) l + 2iV(|p |) <Vi(p)<M-p) 1 ) 

Note that this expression explicitely involves a distribution function, whereas the second equation 
of (H) does not. We will below use ([5]) for a calculation of the physical self-energy.[] Another 
important feature of real-time thermal field-theories is that the effective action has the following 
symmetry || 

rfai,p 2 ] = -r>*,<tf] (7) 

For real, momentum independent vertice (in the present theory with a real scalar field in config- 
uration space), this gives the following relations: 

p(«) _ _r( n _) _ 

1 i\i2-..i n iii2...i n \ ) 

where i = 2 for i = 1 and vice versa. Furthermore one may introduce a functional T[4>] by 

5T[(f)] 5T 



6(f) 5(fi 



(9) 

<Pl=4>,'P2=<P2l<f>} 



where y?2[0] is the solution of the field-equation for (p 2 - For momentum- independent vertice, it is 
easy to see that one has 



ij=l,2 z,j,k=l,2 



and so on. These relations trivially generalize to T\ and we will make use of them below. 

Let us now turn to the calculation of the imaginary part of the self-energy of (at T = 0) 
massless ip -theory. In order to obtain an equation for the scale-dependence of XTI, we plug the 
flow-equations for the two-point functions with 1-1 and 1-2 external legs into the second relation 
in (|5|). We now make a crucial approximation: We neglect the imaginary part of the self-energy on 
the right hand side of all flow-equations (for the present work we also assume 72.fi to be momentum- 
independent, which is however not crucial for the results discussed and could be straightforwardly 
improved). This may be viewed as a " quasiparticle-approximation" . As we will see below, it 
does not influence the leading behaviour for small couplings. Making this approximation however 
means that we will not be resumming imaginary parts in the calculations - a fact important to 
keep in mind. 

In this case, the kernel takes the form HI 



K A (k) = -27i6(k 2 -ml)AS(\k\-A)N(\k \)[ (11 




'In 1, (1 was used for the calculation of the rate. We believe for a number of reasons that will become clear 
below that it is preferrable to use (||). Nevertheless, in principle for an exact solution of the TRG, the outcome 
should be identical. 
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(m\ is the real part of the self-energy) and is purely real. It is thus immediatly clear that all 
contributions to the flow of III are from the imaginary part of the full four-point function. Also, 
due to the specific combination of T^ n and 1^ 12 appearing in (|5|), taking the trace over thermal 
indice in the flow-equation for IT leaves us with the following result: 



dUL A (p) A 3 N{u A ) _ (4 ) 

— = -7-5 e (Po) 2^ ^Aiijk(p-P,QA,-QA) (12) 

OA 47T UJ A i a i, =1 9 



In (|I2D we use the notation 



co A = J A 2 + m\ , Q A = K, |Q A | = A) (13) 



For the quantities appearing on the right hand side of (|12D, we again need the RG-equations 
governing their scale-dependence. Let us first discuss the flow-equation for the coupling appearing 



in (|12D. We use the notation 



1 T A^jk(Pi -P, 9, -q) = iFa\p, -p, q, -<?) (14) 

i,j,k=l,2 

which is suggested by (|TT]|). We now make a further approximation, which is on equal footing with 
the quasiparticle-approximation done above: We neglect the imaginary parts of all couplings on 
the right hand side of the flow-equation for 2Tyy . This should of course not be done in (12), since 



there one would disregard the leading (and in fact only) contribution. Within this approximation 
we find: 

A^fi 4) (p,-p,g,-g) = -^(ri 4) )Y^^ 2 -<)^(|/?+^(l^|-A)iV(|/ |)x 

x Yl m+P + Q) 2 -<Wo+Po + Qo) (15) 

Q=±q 



Note that (|15l) is a completely well-defined expression. Furthermore the vertex appearing here is 
just the combination defined in (|T0|) , and we have taken into account vertice with both 1- and 
2-external legs. This is only possible for the specific sum in (|TJ). Neglecting vertice with 2- legs 
within the present approach would yield ill-defined products of distributions. We have checked 
that the inclusion of trilinear couplings in the present approximations (i.e. considering a theory 
with broken symmetry at low temperatures) is possible along the same line of arguments. Thus 
the results depicted in flT2"|) and (|T5]), together with a flow-equation for the real part of the effective 
action at constant fields, constitute a well-defined nonperturbative method for the calculation of 
thermal rates. This is the first important point in the present work. 

Even though the approach using Wilsonian RG-equations is nonperturbative one wants to 
recover the leading perturbative results in situations where perturbation theory is valid. This is a 
nontrivial problem, since XTI in the theory under study vanishes to one-loop order. 

The loop expansion can be reconstructed from Wilsonian RG-equations iteratively, making use 
of the fact that the right hand side of (0) is down by a factor of h compared to the left hand side. 
In order to obtain a result for ITL to order h 2 , we need the right hand side of (|T^) to order h. The 
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quantities appearing on the right hand side of fll2[) are the thermal mass m\ and the imaginary 
part of the four-point function at finite A. This imaginary part is a pure quantum effect and thus 
is itself of order H, so we need not consider corrections to the mass. 

In order to find the imaginary part of the four-point function (which has to be calculated for the 
specific choice of momenta entering in flT2|) and at nonvanishing external scale A), we turn to the 
corresponding flow-equation (|l|). To order h we may neglect all loop corrections to the quantities 
appearing on the right hand side of (0). Thus to this order the couplings and masses are A- 
independent and the integration of (pT5|) may be performed noting that — 5(\l\ — A) = <9a#(|/| — A). 
One finds 

X f i 4) (p, -p, q, -q) = XTxl^ip, -p, q, -q) + 

+^(r {4 0y^^-- 2 )^(l^+- 2 )^(KI-A)iv(|/ol)x 

xE^(l+J) + Q) 2 - m 2 )e(Z + Po + Q ) (16) 

Q=±q 

Plugging this result into the right hand side of the flow-equation for XTI and performing the 
remaining A-integration should then give the desired imaginary part of the self energy to order 

h 2 . 

The result of a perturbative calculation may be taken e.g. from || where one finds 

2 



s , ^ ( f W ) f d 3 k d 3 q 1 



6 J (2ir) 3 (27r) 3 8u k u q uj r 
x^iSipo+uJk+ujg+uJr) - 6(p -uJk-u q -uj r )] X 

x [l+N k +N q +N r +N k N q +N k N r +N q N r ] + 
+ [5(p +oj k +oj q -oj r ) - 5(po-u k -u q +oj r )] [N r +N k N r +N q N r -N k N q ] + 
+ [6(p +uj k -uj q +u; r ) - 6(p -w k +u q -w r )] [N q + N k N q + N r N q -N k N r ] + 

+ [6{p -u k +LO q +u r ) - 6{p +u k -uj q -uj r )) [N k +N k N q +N k N r -N q N r )\ (17) 

where again uj p = ^p 2 + m 2 and r = k + q '+ p. N p is the Bose-distribution with energy u p . The 
part oc 1 is of course the (T = 0)-imaginary part. In the present method this part enters via the 
boundary conditions and is not calculable within the TRG. The other parts however have to come 
out of our result from (|T2|) and CPL5D - 

An important observation at this point is that for the one-loop contribution in (|16|) one has two 
parts, one being the boundary part at A = oo which is given by the one-loop (T = 0)-result and the 
second one being proportional to and representing the contribution from thermal fluctuations. 
If we simply start with the tree-level effective action at T = - that is set the imaginary part 
of the four point function to for A — > oo - and plug the remaining contribution from (1T0) into 



(J12 ) , we only reproduce the contributions to XII which are bilinear in the distribution functions. 
Dropping all (T = 0)-quantum contributions for the boundary value of the flow of n-point functions 
is however routinely done in applications of the TRG (|jl|, 0, 0, |§, 0). Since the thermal damping- 
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rate vanishes to one-loop, one does not reproduce the leading order perturbative result fll7|) and 
will not be able to give quantitatively reliable results for all temperatures. 

To end the discussion of the perturbative result, we thus need the A — > oo or (T = 0)-value of 
the imaginary part of the four-point function to one loop. This is easily found to be 



irtUp.-P,?,-?) = Ef(r8.)7 



2 r d 3 k 1 

X 



Q=±9 



(27r) 3 UJk^k+pA 



x [8(po+Qo-Uk-Uk+p+o) - 5(p +Q +uJ k +u k+p+ Q)) (18) 

Using this result in (|T6| ) and plugging the resulting expression for TT^ (j>, —p, Qa, —Qa) into the 
flow-equation for T\\{p) given in ( |T2"D one may indeed do the integration with respect to A to find 
for A = the perturbative result ( |17|) up to its (T = O)-part.0 

In order to go beyond the leading order two-loop result, we will take into account the thermal 
renormalization of the real parts of both the self-energy as well as the momentum-independent 
parts of the vertice. We perform a derivative expansion of the full effective action Ta and approx- 
imate it by the effective potential and a standard kinetic term. This approximation is discussed 
in detail in |3|, |7|, || |9] where also the resulting flow-equations may be found. We solve the 
flow-equation for the effective potential as discussed in [0, || |£j and thus resum the momentum- 
independent (thermal) parts of n-point functions with an arbitrary number of external legs - this 
resummation goes considerably beyond the usual daisy- or superdaisy-resummations. It can be 
systematically improved by relaxing the derivative-approximation. 

As the boundary condition for the flow of the effective potential we will use the tree-level 
effective potential of a massless <y9 4 -theory, Ua =00 = a ^- ( t )2 - This amounts to neglecting all (T = 0)- 
loop corrections to the effective potential. Together with the boundary conditions and flow- 
equations for the imaginary parts as discussed above we thus have a closed system of differential 
equations^ which may be solved numerically. 

The resulting imaginary part of the self-energy is connected to the thermal damping-rate j(T) 
through 

= lti A=0 (p = m T ,p = 0) 
2mr 



Perturbatively, using daisy-resummed perturbation theory one finds for a massless theory H, 10, 11 



3/2 

9t=o 



lpeA j "64v/24vrV ^ 



1/2 

1 - ^I^T (20) 



In figure 1 we display the ratio of the result obtained using the flow-equations and (p0|) as a function 
of Qt=o- Indeed the TRG reproduces daisy-resummed perturbation theory for gr=o —> 0, even 

2 In order to obtain the explicit form given in (0), symmetrization in the three- momenta is necessary. The 
calculation is however straightforward. 

3 Note that for the solution of the flow-equation for XT ^ (p, —p, q, —q), the external momentum q is not given by 
(|l3|) with the running values of A and oja- Instead it has to be considered fixed for some independent combination 
(A> A ,). 
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Figure 1: ^(T) / ^ pert (T) as function of log g T =o- 

though the leading perturbative result is two-loop. For finite values of the coupling, the different 
resummations yield different results. Let us again point out that the result of the renormalization 
group-calculations performed here can be understood as a fully resummed result in leading order in 
an expansion in the anomalous dimension and the quasiparticle approximation. It goes beyond the 
resummations used in the literature (it in particular trivially includes the daisy- and superdaisy- 
schemes M). 

Figure 2 shows the dependence of the real- and imaginary parts of the self-energy on the 
external scale A at fixed gr=o- We note the typical behaviour: Small corrections to the (T = 0)- 
values for large A/T due to Boltzmann-suppression of the thermal fluctuations and saturation of 
the values for A <C m^. The limit A — > is completely safe. 

We close this letter by commenting on one of the most important qualitative points made 
in Jl|, namely that a calculation of the damping-rate from the TRG reproduces critical slowing 
down, i.e. the fact that the rate vanishes as the critical temperature of a theory with (at T = 0) 
spontaneously broken ^-symmetry is approached. We believe that even though this conclusion 
was reached on the basis of an incomplete calculation it still holds if one uses the consistent 
approach laid out in the present paper. The reason why this should be the case is easy to see: The 
rate is an on-shell quantity with vanishing external momentum and as the critical temperature 
of the (second-order) phase-transition is approached, the external energy thus vanishes. On the 
other hand, in the framework of Wilsonian RG-equations internal propagators are never massless 
for any A > 0. Close to the phase-transition, the masses and coupling exhibit three-dimensional 
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Figure 2: The thermal mass m\ and the imaginary part of the self energy TYl\(mT) as functions 
of log A/T. 

gr=o = 0.1 was chosen and dimensionful quantities are in units of T. 



scaling behaviour, i.e. the renormalized mass and coupling vanish as 

m T cx(T-T c y ; g T cc (T - T c ) v (21) 

In such a case, the limit uit —>■ should be completely regular and the scaling arguments also 
given in |IJ should still hold: The perturbative result for the imaginary part of the self-energy 
behaves for small as YL(jnT, 0) ~ gT=o^ nm T an d thus one has for the rate 

7 (T)~^lnm T (22) 
ttit 

Consistent resummation using the Wilsonian RG replaces the coupling gT=o with the renormalized 
coupling g T and one obtains with t = T — T c 

7(T^T C ) ocflnt-^0 (23) 

for positive v (in the present model, v is found to be ~ 0.63 |12|, [8|, [!|). 

Within the scheme presented here we now have for the first time a complete approach that 
allows us to study linear-response and static quantities in a nonperturbative manner also in the 
critical regime, while reproducing the perturbative results where they are valid. As shown in 0, 
the TRG can also be formulated for theories involving fermionic degrees of freedom. It should be 
very interesting to apply this nonperturbative method e.g. to questions related to the physics of 
nuclear matter as the chiral phase-transition is approached. 

Acknowledgements: We thank M. Pietroni for helpful discussions on his results. 
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